Strain selection of charge and orbital ordering patterns in half-doped manganites 
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Theoretical and computational results are presented clarifying the role of long-ranged strain in- 
teractions in determining the charge and orbital ordering in colossal magnetoresistance manganites. 
The strain energy contribution is found to be of order 20 — 30 meV/Mn and in particular stabilizes 
the anomalous 'zig-zag chain' order observed in many half-doped manganites. 
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Perovskite manganites of chemical formula 
Lni-^Ak^MnOs (Ln is a lanthanide rare earth such 
as La or Pr; Ak is a divalent alkali such as Sr or Ca) 
present a complex phase diagram associated with the 
interplay of charge, spin, orbital and lattice degrees of 
freedomi Charge and orbitally ordered phases appear in 
wide parameter ranges, for example in Lai-^CaxMnOa 
for x > 0.5, and in Pri_ x Ca x Mn0 3 for 0.3 < x < 0.7. 
The physical mechanism underlying the ordering re- 
mains the subject of debate. Particular attention has 
focussed on the half-doped (a; = 0.5) materials, including 
Lao.sCao.sMnOa, Pro.sCao.sMnOa and Nd .5Sr . 5 MnO3, 
which possess a strongly insulating ground state with 
a particularly interesting ('zig-zag chain') spin, charge 
and orbital ordering pattern shown in the inset to Fig. 
Hi. Stabilizing this state requires 'second neighbor' 
interactions, which as noted by Khomskii and Kugel 2 
and by Ahn and Millis^ may reasonably be expected to 
arise from electron-lattice interactions. However, experi- 
ments have also indicated that long ranged strain plays 
an important role. Insulating ordered phases appearing 
in bulk Lao.5Cao.5Mn03 are suppressed in thin films 4 ; 
it is argued that the difference arises from the lattice 
mismatch with the substrate, which prevents the occur- 
rence of lattice distortions necessary to accommodate 
orbital ordering. Similar conclusions have been drawn 
from recent work on Pr .6Cao.4MnC)3 films subjected 
both to compressive and tensile strain^ Other sets of 
experiments in polycrystalline Lao.sCao.sMnOs show 
that the charge/orbital ordering transition temperature 
is progressively suppressed as grain size is decreased; the 
interpretation again is that boundary effects inhibit the 
formation, in small grains, of the strains imposed by the 
orbital ordering. 6 Perhaps most significantly, polarized 
optical microscopy studies on Bio^Carj.sMnOa single 
crystals^ reveal that as temperature is decreased to just 
below the charge/orbital ordering transition, lenticular 
shaped domains appear and grow slowly with time. This 
behavior is well known in martensitic systems^, where 
it is attributed to the interplay between a tendency of 
the system to favor a state with large strain, and a 



boundary condition which prevents the existence of a 
truly uniform strain. 

Despite the considerable experimental evidence that 
long ranged strain is important, strain physics has 
not received much theoretical attention. The origi- 
nal analyses&ifiiii and many subsequent worksi2iI&I4 fo- 
cused on local interactions, including magnetic, Coulomb 
and Jahn- Teller electron-phonon coupling. In this paper 
we study long ranged strain effects theoretically, and in 
particular show that they play a crucial role in stabilizing 
the particular ordered state characteristic of many of the 
half doped materials. 

The Lni-^Ak^MnOs materials have approximately 
the cubic perovskite crystal structure, in which the Mn 
ions sit on the vertices of a simple cubic lattice. The im- 
portant low-energy electronic degrees of freedom may be 
thought of as Mn e g electrons. There are two e g orbitals 
per Mn. These are degenerate in cubic symmetry. The 
(1 — x) e g electrons are relatively mobile at room tem- 
perature but in many manganite materials, including the 
half-doped compounds of present interest, they become 
strongly localized at low temperatures. The x ^ in- 
sulating states have a charge density which varies from 
Mn to Mn; the pattern of charge density is referred to as 
'charge order'. In many of the half-doped materials, peri- 
odic charge order is observed, with wave vector (jr, tt, 0): 
in other words, in a given plane (which we take to be 
the xy plane) of the ideal cubic perovskite structure the 
charge alternates from site to site in a simple two sublat- 
tice fashion, but the charge ordering simply stacks uni- 
formly in the perpendicular direction. The observation 
of (tt, tt, 0) is surprising: charge ordering typically occurs 
because electrons repel each other, and a fully alternating 
(tt, tt, tt) ordering is typically preferred because it allows 
the greatest distance between charges. 

If charges are localized, then one expects that on a 
given site one linear combination of the two e g orbitals 
will have a higher occupancy than the orthogonal lin- 
ear combination. The pattern of this preferential occu- 
pancy is referred to as 'orbital ordering'. It was first in- 
ferred from the magnetic ordering pattern (discussed be- 
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low) by Wollan and Koehler in the 1950s [T^ using rules 
determined by Goodenough 9 and KanamoriiS, and has 
now been studied in numerous samples by resonant x-ray 
scattering^, which is sensitive both to the orbital occu- 
pancy and to the changes it induces in oxygen positions. 
The orbital ordering occurring in the strongly insulating 
half doped materials may be visualized as 'zig-zag chains' 
on the xy-plane stacked uniformly in the z-direction. The 
'zig-zag' is produced by the alternation of preferred occu- 
pancy of 3a: 2 — r 2 and 3y 2 — r 2 orbitals on the sublattice 
of sites with high charge occupancy (see inset to Fig. [TJi). 
The magnetic ordering pattern corresponds to ferromag- 
netism along the 'zig-zag chains' but with spin direction 
alternating from chain to chain, and follows from the ap- 
plication of the Goodenough/Kanamori rules to the zig- 
zag pattern of statesii. We note, however, that strain 
physics should also be important because the e g orbitals 
have a non-cubic shape. Thus, in the case of 'zig-zag' 
ordering, the occupied orbitals have charge density lying 
mainly in the 'xy' plane, leading to a compression in the 
perpendicular direction. 

We now turn to the theory. A general model would in- 
clude electronic dynamics (hopping and electron-electron 
interaction), electron-lattice coupling, and the energet- 
ics associated with the lattice degrees of freedom. Be- 
cause we wish to focus on the new physics associated 
with strain effects, and because the main application we 
have in mind is to strongly insulating half-doped man- 
ganitcs, for simplicity we adopt a classical, insulating 
model in which the electrons are taken to be localized 
on Mn sites, but our approach may easily be extended 
to an itinerant electron model. In the insulating model 
the e g orbitals of a Mn site can be either empty, or oc- 
cupied by one electron which resides in a linear com- 
bination of the two e g orbitals. We parameterize the 
charge state of site i by a variable hi with hi = for an 
occupied site and hi = 1 for an empty site, so doping 
x = (J2i hi) /N and we focus here oni = 0.5. The or- 
bital state of an occupied site is described by an angle < 
6i < 7r so that \ipi(9i)) — cos8i\d 3z 2_ r 2) + sin9i\d x 2_ y 2). 
(Note that \d^ z 2_ r 2), \d 3x 2_ r 2) l \d 3y 2_ r 2) correspond to 
9i = 0,7r/3,27r/3 respectively while \d z 2_ x 2), |4r 2 -2/ 2 )> 
\d z 2_ y 2), correspond to 9i — 7r/6, 7t/2, 5tv/6.) 

The electron-lattice coupling is taken to have the form 
given in Ref. |13j : a site containing an electron is Jahn- 
Teller coupled to the surrounding oxygen octahedron and 
an empty site induces a breathing distortion of the oc- 
tahedron. In an isolated MnOg octahedron the Jahn- 
Teller effect would lead to an energy gain Ejt = A 2 /2K\ 
with A the Jahn- Teller coupling and K\ the Mn-0 bond- 
stretching force constant (see the inset to Fig. QJ)). Opti- 
cal measurements on LaMnOa suggest K\ w 7.36eV/A 2 
and both band structure calculations and fitting to the 
observed structure suggest A w 1.38eV/A so that Ejt ~ 
0.26eV. We adopt this as the basic energy scale of the 
theory. We write the energy gain due to a breathing dis- 
tortion as Eb = [3Ejt- We expect that the energy costs 
of the two distortions are similar, so /3 ~ 1. 
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FIG. 1: Main panels: ground-state phase diagrams calculated 
by Monte-Carlo annealing on 8 x 8 x 8 lattices, with (a) and 
without (b) the inclusion of the long range energy. The wave 
vector of the charge order (CO) and the pattern 6i of the 
orbital ordering (OO) are indicated. 'Disordered' corresponds 
to phases which show no obvious periodicity or symmetry on 
the lattices studied. Calculations are performed for f3 = 0.8 
for various force constants consistent with bulk modulus Cb = 
0.35-Ki and are presented as functions of the shear C44 and 
Jahn- Teller C* elastic moduli implied by the force constants 
and expressed in units of K\. The inset to Fig. shows 
the 'zig-zag' pattern. Black dots represent the sites with low 
charge occupancy. (b)The 'zig-zag' pattern is absent when 
uniform strain energy is neglected. The inset to Fig. ^> shows 
three of the force constant model couplings. Additional non- 
central force coupling to Mn-Mn-Mn bond angles not shown. 



In the actual perovskite lattice, each O ion is bonded 
strongly to two Mn, so a breathing or Jahn- Teller dis- 
tortion about one site is coupled to distortions at neigh- 
boring sites. To incorporate this physics we follow Refs. 
|13ll7| and write a model whose variables are the devia- 
tions 5i of the Mn ions from their ideal cubic perovskite 
positions and the displacements Ui of the O-ion at z + a/2 
along a direction (a = x,y,z). These variables are cou- 
pled harmonically by near neighbor force constants. The 
inset to Fig. ^> shows the planar force constants. We pa- 
rameterize these force constants by the bulk (Cb), cubic- 
tetragonal (C*) and shear (C44) moduli which they im- 
ply. We then integrate out the variables describing the 
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O positions and obtain an energy of the form 

Etot({hi, 6i,5i}) = £lattice(|oi}) + EjT .breath 

(A detailed derivation will be given in a separate 
communication^) . 

For each electronic configuration ({hi, 9i\) we may now 
obtain the total energy by minimizing Eq. ^ wr th respect 
to the lattice variables We have used a classical Monte 
Carlo simulation process on 4 x 4 x 4 and 8x8x8 Mn 
lattices to find the ground state by annealing the sys- 
tem starting from a random high temperature configu- 
ration. However, an issue arises in the evaluation of the 
energy. To describe this, note that the term Ejt. breath 
may be conveniently written in momentum space as a 
linear coupling between the Fourier transform of the dis- 
placements, 6k, and a vector constructed from the 
Oi and the parameter 0: E JT brcath = e_fcA(fc)£fc. 
The purely lattice term may similarly be written as 
^lattice = \ Y^,k ^-feK(fc)<5fe with K a force-constant ma- 
trix, so that after minimization over the d we have 

E ™ = ~\Y. g - k [A(fc)K- 1 (fc)A(fc)] 4 (2) 
k 

On a finite lattice the sum is over a discrete set of k 
points and the terms with k ^ may straightforwardly 
be evaluated. However, as k — ► we find A(fc) ~ k and 
K(fc) - k 2 so that the lim fe ^ [A(fc)K- 1 (fc)A(fc)] ~ §. 
The limit, and therefore the value of the k — term in 
Eq. [3 depends on the direction from which k approaches 
0. These different limits correspond to different uniform 
strain states of the material (and therefore to different 
changes to shape of the whole lattice), and require spe- 
cial attention. We find that for almost every configura- 
tion of hi,9i there is one particular strain state (direc- 
tion along which the k — > limit is taken) which min- 
imizes E t ot- We use two annealing procedures. In one, 
we use the Etot corresponding to this 'optimum' strain 
state, thereby updating the strain every time the elec- 
tronic configuration is changed and allowing the system 
to minimize its energy with respect to strain, -9 J n the 
other, we neglect the k = term in Eq. ^entirely. One 
could also constrain the k = Fourier components of Si 
to take particular values, thereby forcing the system to 
accommodate to a definite strain (to simulate, e.g. a film 
grown on a lattice- mismatched substrate), but we do not 
explore this physics here. 

We have conducted extensive simulations of Eq. 2] for 
different parameters. The bulk modulus has been esti- 
mated to be 135 GPa (0.35 in units of Ki)2S and we 
consider only choices of Ki consistent with this value. 
We present our results (obtained by varying the Ki in a 
manner consistent with the experimental bulk modulus) 
in terms of the value of the other elastic moduli. Dif- 
ferent values for (3 have been explored, but here we only 
show results for j3 — 0.8. 

Fig. ^shows typical ground-state phase diagrams aris- 
ing from our simulations. The main panel in Fig. 



shows the phase diagram following from optimizing the 
full model over strain. Three states are observed: one 
(CO = (tt,tt, 7r); OO = (7r/4, 37r/4)) corresponds to the 
arrangement of charges and orbitals which implies the 
lowest uniform strain; another, labeled 'Disordered' in 
the figure, corresponds to states which have no obvious 
ordering pattern within the 8x8x8 unit cells we can ac- 
cess numerically, and which presumably either have some 
extremely long-ranged periodicity or are glassy. We note 
that in this parameter regime these 'disordered' states 
are found by the annealing procedure to have energies 
lower than those of any simple periodic structure we have 
devised. 

In the lower right region of the phase diagram in 
Fig. the preferred state corresponds to the one found 
in experiments. Because it involves quite substantial 
anisotropic lattice distortions it is stabilized only in par- 
ticular regions of the phase diagram. The main panel of 
Fig. ^) shows the results of computations for the same 
parameters, but with the uniform strain (k = 0) contribu- 
tion to the energy neglected. We see that the experimen- 
tally observed phase is entirely absent. This conclusion 
remains true in the range 0.6 < (3 < 1.2. Details on the 
role of f3 will be given in a separate communication^ 

The difference in phase diagrams computed with and 
without the strain contribution to the energy suggests 
the importance of long ranged strain effects. To under- 
stand this result in more detail we have evaluated both 
the long range (k — 0) and short range (k ^ 0) con- 
tributions to E to t for many different orbital orders, and 
(tt, tt, 0) and (tt, tt, it) charge order. Typical results are 
presented in Table I. One sees that the differences in 
strain energy between different orbitally ordered phases 
(~ 20meV/Mn) are of the same order as magnetic en- 
ergy differences computed by other workers 14 . Observe 
that the uniform strain energy only depends on the or- 
bital configuration, and is optimized for orbitals which 
correspond to mainly planar charge densities. The 'zig- 
zag' pattern of orbitals is consistent with this result. The 
'uniform clover' pattern (all orbitals 9 = it/2 i.e. d x 2_ y 2) 
has an even lower strain energy, but in the insulating 
model considered here its short range energy is never 
large enough to stabilize this phase. We note, however, 
that electron banding effects are found theoretically to 
favor this phase 21 and that it is observed in metallic 
Lao.sSro.sMnOa. On the contrary, the (ir/4, 37r/4) or- 
bital state is the least strained configuration of all the 
ordered ones, the reason for this being its isotropic charge 
distribution. Note also that the strain term does not dis- 
tinguish between different z-direction orderings: it is the 
short range part of the energy which chooses the stacking 
of charge in the z-direction ((tt, tt, 0) vs (tt,tt,tt)). 

We now comment on possible extensions of our model. 
Magnetic interactions have been neglected, but could be 
easily included. The t2 g localized spins are antifcrro- 
magnetically coupled via superexchange. There is also 
a ferromagnetic superexchange interaction between the 
e g spins which depends on the occupied orbitals in two 
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TABLE I: Energies for different orbital (00) and charge ordering (CO) patterns as indicated, calculated for parameters /3 = 0.8, 
C B = 135 GPa (0.35ifi), C* = 17.35 GPa (0.045/G), and C44 = 30.85 GPa (O.O8.K1) chosen within the range in which the 
'zig-zag' ordering is the ground state (see Fig^i). 





Uniform Strain Energy (meV/Mn) 


Short Range Energy (meV/Mn) 


Total Energy (meV/Mn) 


CO(7r,7r,0), OO 'zig-zag' 


-33.2 


-278.4 


-311.6 


CO(7r,7T,7r), OO 'zig-zag' 


-33.2 


-270.9 


-304.2 


CO(7r,7T,0), 00(7r/2,7r/2) 


-65.0 


-218.0 


-283.0 


CO(7r,7r,7r), 00(7r/4,37r/4) 


-15.9 


-290.8 


-306.7 



neighboring sites. 22 If applied to our 'zig-zag' orbital con- 
figuration, CE-type antiferromagnetism is found. Mag- 
netic energy differences per site between CE-type antifer- 
romagnetism and ferromagnetic ordering are ~ 20meV, 
therefore of the same order than the uniform strain ener- 
gies. Yunoki and co-workers^ have shown how the mag- 
netic interaction may help to stabilize (tt, 7r, 0) ordering 
relative to (tt, tt, tt) ordering; the energetics found in their 
work are sufficiently similar to those found by us that it 
is not clear which mechanism is dominant. Similarly, 
we have adopted a strictly classical, insulating model for 
the electronic degrees of freedom, but our electron-lattice 
coupling and lattice interaction could be easily carried 
over to an itinerant electron model (at least within the 
conventional adiabatic approximation for electron-lattice 
coupling). Use of such a model to explore strain effects on 
a second order weak amplitude ('Charge Density Wave' 
instability) electronic charge and orbital ordering transi- 
tion would be very interesting. Finally, our calculation 



involves only relatively small length scales, and there- 
fore has not included elastic compatibility and other long 
range effects. We hope our results may be useful as input 
parameters for calculations of longer length scale phe- 
nomena. 

In summary, a classical model for charge and orbital or- 
dering has been formulated and analyzed. The important 
new feature is that elastic energies as well as the coopera- 
tive Jahn- Teller coupling are included. Our main finding 
is that energetics associated with long ranged strain are 
at least as important as those from other (e.g. magnetic) 
sources in determining the ground state of the system. 
We find in particular that uniform strain is of vital im- 
portance in stabilizing the 'zig-zag chain' orbital ordering 
pattern observed in many manganite systems. 
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